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The correlator of the square of the Yang-Mills field- strength tensor corresponds to a scalar glue- 
ball, i. e., to a bound-state formed by gluonic ingredients only. It has quantum numbers 
and its mass, as predicted by different theoretical approaches, is expected to lie between 1 and 
2 GeV. Here we restrict our considerations to the Bom level, that is, we consider the correlator 
to zeroth order in the coupling. Gluonic self-interaction is taken into account indirectly by using 
non-perturbative gluon propagators. The employed closed expressions are motivated by lattice 
and Dyson-Schwinger studies. The analytic continuation of the integrals themselves is compli- 
cated by additional obstructive structures like branch cuts and poles that are induced by the inner 
integral in the complex plane of the outer integration variable. We deal with this problem by de- 
forming the outer integration contour accordingly. For different input gluon propagators we find a 
positive glueball spectral density which is required for physical states. Poles are, however, absent 
which is most likely an artifact of working at Bom level. 
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1. Introduction 

As color-carrying states gluons do not appear as asymptotic physical states. They are confined 
to observable color singlet objects by some mechanism, see, e. g., for a short review. In pure 
Yang-Mills theory the only possibility to generate color neutral - and thus observable - states is to 
combine several gluons to form a bound state, a glueball. Experimentally, they are very hard to find 
due to mixing with mesonic states. On the theoretical side, several approaches for gluonic bound 
states are available, see, for instance, for a recent review. 

In the following we calculate the correlator of a glueball candidate and determine its ana- 
lytic properties from which we can extract the spectral density. For a physically observable state 
the spectral density must be positive in order to allow a probabilistic interpretation Q]. Thus, 
the detection of positivity violations indicates that a certain state is expelled form the asymptotic 
state space (and in this sense confined), while the converse does not hold necessarily. Glueballs 
must therefore possess a positive spectral density, although their constituents may not. Indeed, pos- 
itivity violations of gluons are established in the Landau gauge from lattice 0] and functional 
calculations [|, ^. 

As a first approximation we take into account only the Abelian part of the field strength tensor. 
This simplifies the calculations in several respects. For example, renormalization would become 
more complicated but the required machinery is available However, the input we use was 

obtained from full Yang-Mills theory and therefore contains interactions. In the following we will 
use two different fits for the gluon propagator and calculate the glueball correlator numerically. In 



simple cases this is also possible analytically, see, e. g. [|11|, |12p. For future applications a numeric 
procedure is certainly advantageous as it allows, for instance, to use also numerical data as became 
available only recently [p^]. 

The calculation of the correlator boils down to a two dimensional integral. Since we consider 
complex external momenta, an additional subtlety arises: The inner integral leads to non-analytic 
structures, like branch cuts, in the complex plane of the outer integration variable which have to be 
taken into account properly. We do this here by deforming the integration contour, i. e., also the 
radial integration variable becomes complex. This method is already known, see, for example, 
but in the present case the structure of the arising non-analyticities is especially tedious as detailed 
below. 

The two gluon propagator fits we employ here are of the decoupling [jl^, |15|, |16|, ^, |l7|] and scal- 
ing type For the former we use a fit to lattice data motivated by the refined Gribov-Zwanziger 



scenario (RGZ) [|19|]. For the chosen parameter values it possesses two complex conjugate poles. 



Note that other fits, for instance, in [|20[], were suggested as well, but most of them are special cases 
of the used one. For the scaling type propagator we use a fit to the solution of a Dyson-Schwinger 
study It has a branch cut on the negative real axis. 

2. Some Prerequisites 

We consider the correlator of a candidate for a scalar glueball with quantum numbers 

{F\x)F\0)), = {F^,y{x)F^y{x)F^,{0)F^,{OU (2.1) 
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where d is the space-time dimension and F^y{x) is, as part of our approximation, just the Abelian 
part of the Yang-Mills field- strength tensor given by 



We are interested in the momentum space representation of this correlator. 



(2.2) 



(2.3) 



The desired expression, ffd{p ), reads [112] 



(2.4) 



For a transverse gluon propagator we have 



PuPv 



(2.5) 



where only the scalar part ^(p^) enters the expression 1^^. A further complication we have not 
addressed so far is the fact that in 4 Euclidean space time dimensions, the integral as given in 
eq. (|2!4|) diverges like ~ p^. To render the integral finite we employ the BPHZ renormalization, 



i. e., we Taylor subtract the divergent terms: 



^'Ap^)-Mp')-Mo)-p' 



dp'2 



_ 4il 



p=0 



(2.6) 



The odd derivatives vanish because of the anti-symmetry of the angular integral. In order to obtain 
the analytic structure of the scalar glueball correlator, we have to solve eq. ( ^!^ for complex values 
of the square of the external momentum. The spectral density is then accessible by evaluating the 
discontinuity of the branch cut along the negative real axis. For the two-point function A{p^) of a 
given spin zero operator <I>, the spectral density reads 



Pip' 



1 



: lim [A{-p^ -i£)-A{-p^ + i£)] 



(2.7) 



and the spectral representation of the two-point function is 



{2ny 



dT 



T + z' 



(2.8) 



if no poles or cuts except for time-like momenta exist. Tq is the multi-particle threshold. 



Eq. ( p!4j ) holds for arbitrary dimensions. Here we consider only d — 4. The two-dimensional 
case, which of course has a trivial glueball spectrum, served as a test-case for the development of 
the numerics and is presented together with the four-dimensional results in |^. 
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3. The Method 



The algorithm we use here is described in detail in [21], where as an example the analytical 



results from [12] were reproduced. Let us consider the case of the RGZ propagator fit of [|19[], 



C 



;?4 + u^p^ + 1^ ' 



(3.1) 



The fit-parameters are s = 2.508 GeV^ t = 0.72GeV^ u = 0.768 GeV and C = 0.784 [[T|]. In ^] 



the following steps are given in order to evaluate the integral ( |2.4| ) 



• STEP 1: Express ( j2.^ j in hyper-spherical coordinates 



^ .,2^,1, , ,2 [{x + y-2V^^z)y + 2i^V^z-y)^] , (3.2) 
y ~\~ li y ~\~ T 

where x = p^,y = and p-k = ^/x^z. 

• STEP 2: Renormalization 



The integral (3.2) diverges quadratically in x. The renormalized expression is given by ( 



• STEP 3: Analytic continuation 

For the present case this step can be performed either analytically or numerically. For x E C 



the inner integral of eq. ( p^ can produce an integrable singularity together with the rest of 
the integrand. When z runs through its integration interval [—1, 1], it picks up a whole line of 
these singular points resulting in a branch cut in the complex plane of the radial integration 
variable y. Thus the contour of the radial integral has to be deformed in order to avoid the 
cut. For eq. (|^ we find two branch cuts as well as a pair of complex conjugate poles. The 



branch cuts, parametrized by z, in the y-plane can be determined analytically by finding the 



zeroes of the integrand of eq. ( p.2| ) for a given x G C. We compared these results with a 
numerical integration. For x = — 2 + 2/ both are shown in Fig. |l|. 

It is clearly visible in Fig. |l] that the deformation of the contour of the y-integration, required 
to connect y = toy = where is a UV cutoff, can be quite tricky. In general the open 
piece between the branch cuts always points in the direction of Arg(x). Thus, if x is on the 
positive real axis, the integration is straightforward since 3; can be kept real as well. Now 
let us consider a complex x = (r, 0) by keeping r fixed while < < 7r/2. There are no 
poles in the first quadrant, and the opening of the branch cuts always point in the direction 
of Arg(x), thus the contour can be deformed continuously in that case. The same is true 
for the fourth quadrant. However, the complex conjugate poles of the integrand located in 
the quadrants II and III require more care. Obviously the contour cannot be deformed as 
easily for Arg{pui) > Arg(x) > Arg(/7//), where pu and puj are the pole locations in the 
second and third quadrants, respectively. For some values of x the branch cut end points are 
narrowing down the area for a possible contour, see Fig. |[ When an endpoint of a branch 
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Figure 1: Left: Analytic results for the branch cuts and poles in the complex _y-plane for x = — 2 + 2i. Right: 
Numerical verification of the analytic result. 
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Figure 2: Left: The region for possible contours narrows down, x = —2. Right: x = — 1 . 18 + 2.7/. The same 
situation occurs for x = — 1 . 18 — 2.7/, with the cuts flipped around the real axis. 



cut coincides with one of the poles, the contour cannot be deformed continuously and a non- 
analyticity arises in the integral. In we confirmed that for all points where this happens a 
branch point is also predicted from the Cutkosky rules [E2[]. 



There are two further steps which we omit here as they are purely technical. What is relevant 
here is that the complex conjugate poles together with the two branch cuts severely restrict the 
possibilities for the contour deformation. It is hard to obtain stable results for complex values of 
X when the argument of x coincides with the argument of one of the pole locations. As discussed 
in the next section, we find three branch cuts for the RGZ case, one along the negative real axis, 
and two along the directions Arg(p//) and Arg(p///). The numerical determination of the branch 
points in this case is very troublesome, because for x- values close to the cuts in the x-plane the 
contour necessarily always comes very close to the cuts in the y-plane what leads to numerical 
artifacts. Even though the scaling propagator of [|[] has a branch cut and the integrand induces two 
more cuts in the y-plane, the absence of poles allows a continuous contour deformation to values 
very close to the negative real axis. The results for the scaling propagator are thus not plagued by 
numerical issues. 
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-5^5 -5^5 



Figure 3: Left: The imaginary part of the scalar glueball correlator with RGZ gluons as input. Right: The 
real part of the correlator. 

4. Results 
4.1 Decoupling 

In the previous section we already discussed several aspects of the RGZ propagator as gluonic 
input. Most importantly, we confirmed the location of the branch points known from the Cutkosky 
rules. Fig. || shows the imaginary and real parts of the correlator. The three branch cuts are clearly 
visible. The two 'unphysical' ones open very slowly. The extracted discontinuity of the 'physical' 
branch cut is depicted in fig. Q It becomes negative for small values of —p^ and rises earlier than 
expect from the Cutkosky analysis. From investigating the complex plane of the radial integration 
variable we know that these phenomena are numerical artifacts which we expect to vanish if the 
contour deformation is better tuned; see ^ for a more detailed discussion. Thus we conclude that 
the spectral density is positive. 




-p' [Gev'] -p'[GeV'] 
Figure 4: The discontinuity of the physical branch cuts. Left: Decoupling gluons. Right: Scaling gluons. 
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Figure 5: Left: The imaginary part of the scalar glueball correlator with scaling gluons as input. Right: The 
real part of the correlator. 

4.2 Scaling 

The IR part of the scaling gluon fit is given by [|8|], 

J / 2 

^('^) = -J^y^) . <"■') 

with K = 0.595353. The exponent leads to a branch cut of the propagator for complex momenta. We 
neglected the UV part of the propagator fit, which involves a logarithm, as we are only interested 
in IR relevant parts of the propagator. The other parameters are w = 2.5 and A = 0.51 GeV. 

Fig. H shows the imaginary and the real parts of the correlator. Since there are no non- 
analyticities besides the branch cut on the negative real axis, a Kallen-Lehmann representation 
is possible. The corresponding positive spectral density is depicted in fig. We also observe that 
in this case the evaluation of the correlator in the complex plane is not plagued by numerical ar- 
tifacts. Strictly speaking the Cutkosky analysis is in this case not applicable, since the employed 
propagator does not have the required form. However, a naive application leads to a threshold in 
precise agreement with our numeric result. 



5. Summary 

In this work we studied the analytic properties of a scalar glueball correlator at Bom-level. The 
self-interaction of gluons entered via using non-perturbative gluon propagator fits. These exhibit 
positivity violations and describe thus confined gluons. The resulting glueball correlators have a 
branch cut for time-like momenta and no poles. The extracted spectral densities are positive as 
required for a physical state. For the fit of the decoupling propagator we also find two unphysical 
cuts which are due to the analytic structure of the fit. Possible continuations include the addition 
of higher order terms and the use of numerical results for the propagator in the complex plane. 
The employed techniques for the evaluation of the integrals may be useful for other studies, where 
complex momenta are involved, as well. 
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